Quasiparticle entropy in the high-field superconducting phase of CeCoIn 5 
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The heavy-fermion superconductor CeCoIns displays an additional transition within its super- 
conducting (SC) state, whose nature is characterized by high-precision studies of the isothermal 
field dependence of the entropy, derived from combined specific heat and magnetocaloric effect mea- 
surements at temperatures T > 100 mK and fields H < 12 T aligned along different directions. 
For any of these conditions, we do not observe an additional entropy contribution upon tuning at 
constant temperature by magnetic field from the homogeneous SC into the presumed Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) SC state. By contrast, for H || [100] a reduction of entropy was found 
which quantitatively agrees with the expectation for spin-density-wave (SDW) order without FFLO 
superconductivity. Our data exclude the formation of a FFLO state in CeCoIn5 for out-of-plane 
field directions, where no SDW order exists. 



Superconductivity can be affected by an imbalance of 
density of states (DOS) for spin-up and spin-down elec- 
trons introduced by Zeeman splitting under magnetic 
field, leading to a finite total momentum Cooper pair- 
ing, in contrast to the usual BCS pairing. Since Fulde, 
Ferrell, Larkin and Ovchinnikov [lj, |2[ (FFLO) have pre- 
dicted such a superconducting (SC) state more than forty 
years ago, numerous attempts for its experimental re- 
alization have been made in various systems, including 
heavy fermion (HF) superconductors [3(, organic super- 
conductors [j and cold atoms on optical lattices [5[ . HF 
superconductors are promising due to their largely en- 
hanced Pauli paramagnetic susceptibility. 

The tetragonal HF compound CeCoIn 5 undergoes a SC 
transition at 2.3 K at ambient pressure and zero magnetic 
field [6(. Power-law behavior in the low-temperature 
specific heat and thermal conductivity suggests an un- 
conventional SC state [7J, while the extraordinary large 
electron mean-free path of ~4000A [8| sets SC in the 
clean limit. The Maki parameter, which characterizes 
the relative strength of Pauli to orbital limiting effects 
in magnetic field, ranges between 3.5 (H | [001]) and 
4.5 (H | [100]) [9j, i.e. strongly exceeds the minimal 
value of 1.8, required for the FFLO state [10]. The 
discovery of an additional phase in the high-field and 
low- temperature (HFLT) corner of the SC phase dia- 
gram [ll|, [l2j was therefore proposed to be the real- 
ization of the long sought-after FFLO state and pro- 
moted numerous experimental and theoretical studies. 
Specific heat shows an anomaly across the transition 
between the low-field SC and the HFLT SC phases at 
magnetic fields above -Hhflt = 10.3 T along the tetrag- 
onal basal plane (H c2 = 11.7 T) [HI, El- This HFLT 
SC phase was further confirmed by several measurement 
techniques [l3l4l9j . However, the formation of a FFLO 
phase in this material is controversial, since the sub- 
sequent NMR and neutron scattering studies found an 
incommensurate small-moment antiferromagnetic (AF) 
order in the HFLT SC state, which is likely of spin- 



density- wave (SDW) type [2fJ, |2l|. One of the most 
peculiar properties of this phase is that the AF order 
does not extend into the normal state and exists only in 
the SC state |2l|, suggesting some additional stabiliza- 
tion of AF order by the SC state. Some theories pro- 
posed mechanisms for stabilizing AF order due to strong 
Pauli-limiting and a nodal SC gap structure, without a 
FFLO state 22j, [23J, while in another theory, a coex- 
isting FFLO state is necessary for the formation of AF 
order [24| . A recent In-NMR study for fields along [100] 
has suggested the formation of a pure FFLO phase lead- 
ing to an anomalous line broadening already at fields be- 
tween 9.2 T and i^HFLT = 10.3 T, while the coexistence 
of AF order and FFLO superconductivity is claimed at 
-Hhflt < H < H c2 [25(. A spatially uniform coexis- 
tence of AF order and FFLO nodal planes has also been 
suggested from the most recent NMR study [26[ . The or- 
dered moment associated with the AF order disappears 
when the field is rotated by more than 17° out of the 
basal plane [271 ] . Thus, any remaining anomalies, in par- 
ticular for H || [001], could not be related to AF order. 
The Maki parameter for this field direction is still twice 
as large as the required value for the formation of FFLO 
state, but much less studies on the SC state at high fields 
for H || [001] have been reported. NMR experiments [19J] 
suggest a FFLO state with a rather temperature inde- 
pendent phase boundary around 4.7 T, i.e. very close 
to -Hc2=5.0 T, in contrast to the strongly temperature 
dependent 7?hflt transition along the basal plane direc- 
tion 



11|. 



In the FFLO state, the SC gap function, A(r), is spa- 
tially modulated with a wave length, 2w/q, and param- 
agnetic quasiparticles appear periodically at nodal po- 
sitions (A=0). Because of the additional quasiparticles 
in the FFLO state above the critical field, i^FFLO, the 
isothermal entropy as a function of field shows an addi- 
tional convex t/H — -ffFFLO-like contribution, leading to 
a steep increase of S{H,T = const) [28|, [29(. Any mag- 
netic ordering, by contrast, will have a negative isother- 



mal entropy contribution related to the reduction of de- 
grees of freedom. Therefore, the observation of a steep in- 
crease would be a "smoking gun" proof of FFLO SC. Ex- 
perimentally, we can very precisely determine the isother- 
mal field dependence of the entropy by measuring the 
magnetic Griineisen ratio Tji = ^/T(dT/dH)s, i.e., the 
magnetocaloric effect under perfect adiabatic conditions, 
and the specific heat C using the thermodynamic relation 



dS 

dH 



C dT 
T dH 



(1) 



High quality single crystals were grown by the self-flux 
method. The specific heat and magnetic Griineisen ratio 
were measured with very high resolution in a dilution re- 
frigerator with a SC magnet equipped with an additional 
modulation coil by utilizing heat-pulse and alternating 
field techniques, respectively [30[. Using (1), we could re- 
solve entropy changes as small as 2 x 10~ 5 J/mol-K within 
the SC state of CeCoIns, corresponding to 3.5 ppm of 
i?log2. 

We first concentrate on measurements for H || [100], 
shown in Figure 1. The overall convex shape of the field- 
dependence of the heat capacity at 0.2 K results from the 
strong Pauli limiting effect of the SC 3l| . Pronounced 



peaks in the heat-capacity and magnetic Griineisen ra- 
tio indicate the first-order transition to the normal state 
at H C 2 — 11.5 T. The additional second-order transition 
between the regular and HFLT-SC states at H-rf leads 
to broadened discontinuities in the two quantities. The 
temperature dependence of Hrf (cf. Inset Fig. lb) is 
in perfect agreement with previous results |ll| . Using 
(1), we determine the field-dependence of the entropy at 
0.2 K (cf. Fig. lb). At H H f=10.4 T, the field deriva- 
tive (dS/dH)T displays a broadened downwards discon- 
tinuity indicating a negative contribution to the entropy. 
The overall increase of entropy with field is naturally ex- 
plained by the increasing number of paramagnetic quasi- 
particles. We do not resolve a phase transition at 9.2 T, 
i.e. the field at which a drastic broadening of the In-NMR 
spectra has suggested "exotic superconductivity" [25j. 
This excludes the formation of a FFLO state in this part 
of the phase diagram. 

Since previous neutron diffraction measurements have 
proven that the AF state disappears at out-of-plane an- 
gles larger than 17° 27J, in accordance with recent the- 
oretical work [32j , it is of particular interest to compare 
the measurements for H | [100] with respective mea- 
surements in tilted field, cf. Figure 2a. Indeed the 
thermodynamic signatures of the HFLT transition have 
disappeared at a tilting angle of 18° for which mono- 
tonic behavior up to H C 2 is found in the various prop- 
erties [l6(. Using the field dependence of (dS/dH)x, we 
extract with high precision the isothermal field depen- 
dence of the entropy (inset of Fig. 2b). At H < Hhf 
and H > H C 2 similar behavior is found for both field 




FIG. 1: (Color online) (a) Magnetic field dependence of elec- 
tronic specific heat divided by temperature (black solid cir- 
cles, left axis) and magnetic Griineisen ratio (red open circles, 
right axis) of CeCoIn 5 at 0.2 K for H || [100]. (b) Calculated 
field-derivative of the entropy (red open circles, right axis) and 
integrated entropy increment (black solid circles, left axis). 
Inset displays field-dependence of magnetic Griineisen ratio 
at various different temperatures. Arrows indicate the tran- 
sition to the high-field SC phase, Hhf- 



orientations, while in the intermediate field regime a de- 
pression is visible for H || [100]. This negative entropy 
increment of 8 mJ/mol-K at 0.2 K (cf. two-sided ar- 
row in the inset) must be related to the AF ordering 
for H || [100], which disappears at 18°. Using the value 
of the specific heat coefficient at the same temperature, 
Coi/T = 0.67 J/mol-K 2 , we estimate a 6% reduction of 
the DOS at the Fermi energy. Chromium, for exam- 
ple as comparison, loses a similar fraction, ~4% of DOS 
due to the gapping of the Fermi surface at the SDW 
transition 



331 ] . The fraction of truncated Fermi surface 



area, t, can be estimated from the SDW modulation vec- 
tor, Q=(0.56,0.56,0.5) [Hj], and the SDW gap, A SDW , as 
t=(2pAsDwm*) / (\Q\ n 2 ) where p is the number of trun- 
cated faces and to* the effective mass [33(. The gap is 
approximated by the BCS formula, AgDw=3.5fcBTsDW, 
using the phase transition temperature of 0.27 K ll| . For 
a tetragonal system p — 4 and to* can be estimated from 
the Sommerfeld coefficient, 7=0.67J/mol-K 2 , as m* — 
(3h 7)/(fc 2 ! V m £;F)=230mo for a spherical Fermi surface 
and to* = (3?i 2 7)/(A:|a 2 iVA) = 1000TOo for a cylindrical 
Fermi surface. Here V m is molar volume, fep = Q/2, 
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FIG. 2: (a) Electronic specific heat divided by temperature 
of CeCoIns at 0.2K for H || [100] (solid circles) and 18° tilted 
from [100] towards [001] (open circles) versus magnetic field 
normalized by upper critical field H c2 = 11.7 T (H || [100]) 
and 9.0 T (H tilted by 18° towards [001]), respectively. In- 
set shows the respective magnetic Griineisen ratio data, (b) 
Calculated field-derivative of the entropy. Inset shows respec- 
tive entropy increments for H || [100] (solid circles) and 18° 
(open circles). Solid arrows indicate the high-field transition 
for H || [100]. Two-sided open arrow in inset indicates the 
entropy difference of 8mJ/mol-K. 



a = 4.62A and N\ denotes Avogadro's constant. Then, t 
equals 0.015 and 0.065 for spherical and cylindrical Fermi 
surfaces, respectively. In the quasi-2D case relevant for 
CeCoIns and intermediate value is expected, in good 
agreement with the experimentally observed t — 0.06. 

We also estimate the increase of the entropy for a pos- 
sible FFLO state, which arises due to a spatial mod- 
ulation of the SC gap along the direction of the ap- 
plied field. The most energetically stable wavelength 
of the gap modulation in the FFLO state was theoret- 
ically determined as a function of magnetic field [34j . A 
long wave length right above .Hfflo decreases sharply 
as the field is increased and reaches the smallest value 
of 17 Rq at a magnetic field close to H c2 , where 



Rq=£Jvvf) I '(2TrkBT c ) and vp denotes the Fermi veloc- 
ity J3l|. The reported in-plane v F = 7000 m/s [HIH, 
yields i? =38A and the shortest modulation wavelength, 
17i? =650A. Using the SC gap A =4.92K [7], we can cal- 
culate the spatial dependence of the quasiparticle occupa- 
tion f{z)~[l+eyi\){Ey(z)/k'QT)\~ 1 along the z-direction 
(parallel to the field). Here E k (z) = [(e k - E F ) 2 + 
(Ao sin(27T2;/17i?o)) 2 ] 1 ^ 2 and ek denotes the kinetic quasi- 
particle energy. Note, that at the nodal positions / 
equals the normal Fermi Dirac function, giving rise to 
an enhancement of the entropy compared to the SC 
state. The position dependent entropy is derived from 
S(z) = -2fc B S,[(l - /(*)) ln(l - /(*)) + f(z) ln(/(z))]. 
Its z-average at 0.2 K amounts to 15% of the entropy 
difference between H HF and H c2 . Using S(H c2 ) — 
5(J5 r HF)=25mJ/mol , K, this 15% increment corresponds 
to 4 mJ/mol-K. Consequently, within a SDW-FFLO co- 
existence scenario, this increase needs to be overcompen- 
sated by an entropy reduction due to the SDW ordering 
of 12 mJ/mol-K at 0.2 K, in order to match the experi- 
mental data. This would correspond to t ps 0.09, which is 
well beyond the estimated upper limit for the fraction of 
DOS which could be gapped due to the SDW formation. 
We therefore conclude that the data indicate a SDW or- 
dering without FFLO state. 

We now turn our attention to measurements in mag- 
netic fields applied 90° off the basal plane, i.e. paral- 
lel [001]. Figure 3 shows the isothermal field depen- 
dences of C e j/T, Th and the corresponding (dS/dH)? 
and AS at 0.2 K. The overall concave shape of C e i/T 
is similar as found for H || [100]. A sharp peak in the 
field-derivative of the entropy at H c2 — 4.9 T is char- 
acteristic of the first order SC transition, caused by the 
strong Pauli limiting [9j. Whereas the isothermal field 
derivative of the entropy (dS/dH)x has shown a step- 
like decrease at H | [100] compatible with a second- 
order phase transition, for H || [001] it displays only a 
change in slope. Recent isothermal magnetization mea- 
surements have found related behaviors for the magnetic 
susceptibility dM/dH along the two field directions [37]. 
Changes in slope of second-order derivatives of the free 
energy, such as (dS/dH)? correspond to discontinuities 
in third derivatives of free energy only, i.e. too weak for 
a second-order phase transition. With increasing tem- 
peratures, the observed kink broadens significantly and 
is shifted towards lower field values, as shown by the red 
squares in the inset of Fig. 3(a). (dS/dH)x for the field 
applied 70° tilted away from [100] towards [001] is also 
plotted in the inset of Fig. 3(b), to show the evolution of 
the kink with the field angle. It is still visible, however 
broadened. Presumably, this feature broadens further 
with tilting angle and smoothly changes to a featureless 
curve, similar to that for the field 18° off the basal plane. 

The origin of this " kink-signature" could not be related 
to the previously observed AF ordering, which has shown 
to disappear at angles larger than 17° 27). The forma- 



tion of a FFLO state, suggested by earlier NMR experi- 
ments [ljj, could be excluded as well, since, as explained 
above, a sharp increase of (dS/dH)? upon increasing the 
field across -ffFFLO due- to a sudden increase of the quasi- 
particle entropy would have been expected. Another pos- 
sibility would be that the anomaly is related to a change 
of the vortex lattice. Small angle neutron scattering has 
previously revealed a first-order rhombic to hexagonal 
transition of the vortex lattice near 4.4 T [38[. These ex- 
periments have detected also a second-order square to 
rhombic transition around 3.3 T, which our measure- 
ments do not detect. Observation of vortex lattice phase 
transitions by bulk measurements is very rare, except for 
solid-liquid transitions in cuprates, exhibiting entropic 
anomalies 39|, |40j. 

Previous Hall effect and thermal expansion measure- 
ments in the normal state of CeCoIn5 for H \\ c have 
found crossover lines T cr (H) and Tfl(H) separating non- 
Fermi liquid from Fermi liquid behavior, which display a 
linear temperature vs field relation and extrapolate to 
a critical field of roughly 4 T in the zero-temperature 
limit [411 . |42j . This suggests a quantum critical point 
(QCP), hidden by the SC phase and possibly related 
to the suppression of AF order emerging under negative 



pressure or Cd-doping [42, |43j. Generically, the mag- 



netic Griineisen parameter and thereby (dS/dH)x are 
expected to display a characteristic sign change within 
the quantum critical regime, due to the accumulation 
of entropy 44]. For CeCoIn 5 , the field dependence of 
the entropy in the SC state is dominated by the vortex- 
lattice quasiparticle contribution. Assuming a smooth 
evolution of the latter, as sketched by the dashed line in 
the inset of Fig. 3b, the observed kink-like anomaly may 
result from an additional quantum critical contribution 
which changes sign. Note, that quantum criticality at 
finite temperatures could result from a nearby QCP in 
multi-parameter space, e.g. being located at a finite Cd- 
doping [43j. Finally, we note that the observed anomaly 
further broadens upon tilting the field from [001] towards 
[100] and disappears at 18° from the in-plane field orien- 
tation. 

In conclusion, our high-precision measurements of the 
isothermal field dependence of the entropy in the SC state 
of CeCoIn.5 do not find an additional component, which 
would indicate nodal quasiparticles in a FFLO SC state. 
By contrast, for isothermal measurements at H | [100] a 
clear reduction of the entropy by 8mJ/mol-K at 0.2 K is 
found at a second-order HFLT transition at 10.4 T. This 
transition coincides with the previously detected incom- 

m si 



mensurate AF order 



The observed reduction of 



DOS is in perfect agreement with the expectation for a 
SDW formation without additional FFLO state. Upon 
tilting the field direction by 18° towards the [001] di- 
rection, the HFLT transition has completely disappeared 
and a FFLO state could be excluded within the exper- 
imental resolution, which is three orders of magnitude 




FIG. 3: (Color online) (a) Electronic specific heat divided 
by temperature (black solid circles, left axis) and magnetic 
Griineisen ratio (red open circles, right axis) of CeCoIns as 
a function of field for H || [001] at 0.2 K. Inset shows the su- 
perconducting phase diagram for H || [001]. Black circles and 
red squares mark positions of upper critical field and "kink 
anomaly" (cf. black arrows in main panels), respectively, (b) 
Calculated field-derivative of the entropy, (dS/dH)r, (red 
open circles, right axis) and integrated entropy increment 
(black solid circles, left axis). Inset shows {dS/dH)r for the 
field applied along [001] and different angles tilted away from 
[100] towards [001]. Dashed line see text. 



better than the estimated increase of entropy due a FFLO 
state. A FFLO state could also be excluded for H || [001], 
where a broadened kink anomaly is found, which likely is 
related to a nearby QCP. Finally, we note that a similar 
study of the isothermal field dependence of the entropy 
could provide a conclusive test for the existence of FFLO 
SC states in other candidate materials such as organic 
superconductors [4J, |45|, |46 1 . 
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